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Abstract 

The  computational  fluid  dynamics  (CFD)  study  of  the  mass  and  heat  transfer  in  the  phosphotungstic  acid  electrolyte-based  fuel  cell 
(PWAFC)  was  performed  in  order  to  elucidate  the  cause  of  irregular  fuel  cell  work  failures  at  high  current  densities,  which  manifested 
through  partial  melting  of  the  matrix  and  electrodes  and  crossover  effect.  It  was  established  that  the  reason  for  these  difficulties  is  not  the 
damage  (cracks)  of  the  catalyst  layer  in  the  gas  diffusion  electrode,  but  the  bad  construction  of  the  reactant  gas  and  liquid  electrolyte 
distributors  on  cathode  and  anode  side  of  the  PWAFC  and  the  lack  of  cooling  system.  The  CFD  study  pointed  to  several  practical 
conclusions  how  to  improve  the  PWAFC  construction.  ©  2001  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

A  low-temperature  hydrogen/oxygen  fuel  cell  has  been 
developed  recently,  based  on  the  concentrated  aqueous 
solution  of  phosphotungstic  acid  (PWA)  as  proton  conduct¬ 
ing  electrolyte  [1,2].  Previous  attempts  to  use  PWA  as  a  solid 
proton  conducting  electrolyte  in  crystalline  form  had  failed 
in  spite  of  its  high  proton  conductivity  (0.18  S  cm  ')  at 
room  temperature  and  high  relative  humidity  of  surrounding 
atmosphere  [3-6].  The  solid  PWA  electrolyte  must  be  kept 
in  its  most  hydrated  crystalline  form  (H3PW1204o-29-30 
FLO),  which  has  high  enough  proton  conductivity 
(0.17  S  cm-1)  to  be  applied  in  fuel  cells  [3].  This  crystalline 
form  exists  only  in  the  narrow  interval  of  the  high  relative 
humidity.  If  the  humidity  drops  below  80%,  the  conductivity 
of  solid  PWA  drops  by  nearly  an  order  of  magnitude.  On  the 
other  side,  if  the  humidity  approaches  100%,  the  phase  of 
PWA  changes  abruptly  from  solid,  crystalline  electrolyte  to 
the  concentrated  water  solution  of  PWA,  which  has  high 
proton  conductivity  (over  0.3  S  cm  1  at  295  K)  [7].  More¬ 
over,  the  change  of  the  PWA  electrolyte  phase  occurs  non- 
uniformly  within  the  solid  electrolyte  in  the  working  fuel 
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cell  while  it  depends  on  the  current  density  distribution  on 
the  electrodes.  The  latter  depends  in  its  turn  on  the  gas  flow 
and  distribution  within  the  porous  catalyst  layer  of  the  gas 
diffusion  electrodes.  The  water  management  of  the  fuel  cell 
with  solid  PWA  electrolyte  thus  becomes  a  very  demanding 
engineering  task. 

Since  the  idea  to  use  PWA  as  the  proton  conducting 
electrolyte  was  originally  founded  on  the  ability  of  this 
electrolyte  to  modify  the  Pt  electrode  surface  in  such  a 
manner  that  it  becomes  more  resistant  to  transition  metal 
ions  poisoning  [8]  and  to  the  CO  poisoning  [9]  and  that  it 
enhances  the  kinetics  of  HER  [10],  we  decided  to  continue 
the  development  of  this  new  type  of  fuel  cell  even  by  using 
the  PWA  as  liquid  electrolyte  in  the  form  of  concentrated 
aqueous  solution.  However,  it  was  soon  realized  that  the 
utilization  of  a  traditional  cell  with  the  electrolyte  flow 
parallel  to  the  electrodes  was  not  suitable  for  this  type  of 
electrolyte  [7].  The  high  density  and  high  dynamic  viscosity 
of  the  concentrated  PWA  solution  (60%  PWA  solution  has  a 
density  of  about  2  g  cm-3  and  dynamic  viscosity  of  about 
1 .3  cP  at  room  temperature)  hinder  its  flow  inside  the  porous 
matrix,  creating  high  difference  in  PWA  concentration 
between  the  fuel  cell  inlet  and  outlet  points.  Therefore, 
the  through-flow  configuration  of  the  cell  was  developed 
with  the  flow  of  concentrated  PWA  solution  in  the  direction 
from  cathode  to  anode  [2] .  With  this  configuration  we  have 
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successfully  resolved  also  the  problem  of  mixed  potential 
created  by  electron  conductivity  of  PWA  via  redox  reaction 
[11].  Of  course,  the  preparation  of  the  gas  diffusion  electro¬ 
des  had  to  be  optimized  for  the  use  with  the  concentrated 
aqueous  solution  of  PWA  in  order  to  create  a  proper  three- 
phase  interface  in  which  the  electrochemical  reactions  pro¬ 
ceed  [12].  The  results  obtained  with  this  cell  were  encoura¬ 
ging.  The  40  cm2  cell  with  Pt  loading  of  0.5  mg  cm  2  in  gas 
diffusion  electrodes  working  at  room  temperature  and  nor¬ 
mal  pressure  was  able  to  produce  a  maximum  output  power 
density  of  about  0.7  W  cm-2  at  current  density  1.5  A  cm~2 
[1].  This  performance  was  comparable  to  the  laboratory 
performance  of  proton  exchange  membrane  fuel  cell 
(PEMFC)  with  Dow®  membrane,  having  0.4  mg  cnT2  of 
Pt  and  working  at  323  K  and  normal  pressure  [13], 

Further  development  of  the  proof-of-the-concept  fuel  cell, 
based  on  concentrated  liquid  electrolyte,  has  manifested 
some  serious  drawbacks  in  the  fuel  cell  operation:  the  loss 
of  electrolyte  through  the  outlet  tubing  of  reactant  gases  and 
the  excessive  consumption  of  reactant  gases  with  respect  to 
their  stoichiometric  consumption  in  electrochemical  reac¬ 
tion.  Further  post-mortem  visual  analysis  showed  zones  in 
which  the  layers,  electrodes  and  matrix,  were  dried  and 
melted  together  due  to  hydrogen  cross-over  and  hot  spots 
formation  on  the  surface  of  cathode  and  other  zones  where 
the  electrodes  and  matrix  appeared  flooded  with  acid  elec¬ 
trolyte.  By  these  observations,  it  became  evident  that  there 
were  parts  of  electrodes  and  matrix  assembly  that  operated 
under  different  fluid  dynamic  conditions.  However,  it  was 
not  possible  to  ascribe  this  damage  of  PWAFC  components 
unambiguously  to  one  reason.  They  could  be  a  consequence 
of  inappropriate  procedure  for  fuel  cell  components  assem¬ 
bling  during  which  the  small  cracks  in  the  catalyst  layer  of 
electrodes  could  be  created  and  this  would  lead  to  the 
hydrogen  crossover.  But,  they  could  also  be  created  by 
inappropriate  design  of  the  reactant  gases  and  liquid  elec¬ 
trolyte  flow  patterns. 

The  complex  problems  of  inappropriate  fluid  flow  and 
distribution  within  the  fuel  cell  can  be  resolved  by  using 
computational  fluid  dynamics  (CFD)  approach.  This 
approach  has  an  additional  advantage:  the  detailed  and 
comprehensive  analysis  can  be  applied  to  different  fuel  cell 
geometry  and  can  be  used  simultaneously  also  for  the 
resolution  of  heat  transfer  aspects  of  design.  There  are 
several  examples  of  numerical  simulation  of  fluid  flow 
and  heat  transfer  in  fuel  cells  described  in  the  literature 
[14-18].  Several  two-dimensional  models  of  fluid  flow  in 
PEMFC  have  been  developed  recently  in  order  to  elucidate 
the  influence  of  the  geometry  of  gas  distributors  and  of 
porosity  and  thickness  of  diffusion  layer  of  the  gas  diffusion 
electrode  (mainly  cathode)  on  the  current  distribution  over 
the  surface  of  the  gas  diffusion  electrode  [17-20].  With  the 
advance  of  the  PC  computer  performance  the  rather  sophis¬ 
ticated  software  packages  can  be  used  to  resolve  engineering 
problems  which  would  otherwise  take  to  much  time  and 
efforts  to  be  resolved  experimentally  [21]. 


In  order  to  analyze  the  observed  anomalies  in  the  PWAFC 
performance,  we  decided  to  simulate  the  fluid  flow  and  heat 
transfer  in  a  complete  1:1  physical  model  of  the  PWAFC 
using  CFD  program  PHOENICS,  Version  2.2  for  PC  [22]. 
The  results  of  the  three-dimensional  numerical  simulation 
were  validated  through  gas  flow  experiments  performed  in 
the  cathode  side  half-cell  configuration.  The  detailed  ana¬ 
lysis  of  the  reactant  gases  flow  and  of  the  electrolyte  flow  as 
well  as  the  analysis  of  heat  transfer  within  the  PWAFC  are 
presented  in  this  study.  The  conclusions  were  used  to  design 
PWAFC  with  better  geometry. 


2.  Description  of  the  PWAFC  physical  model 

2.1.  PWAFC  geometry 

This  type  of  fuel  cell  configured  as  a  through-flow  porous 
catalytic  plate  reactor  [2]  is  composed  of  several  square  or 
rectangular  thin  plates  made  of  different  electron  conducting 
and  insulating  materials.  This  thin  plate  technology  was 
used  deliberately  to  simplify  the  production  of  all  elements 
and  the  assembling  of  the  fuel  cell  and  fuel  cell  stack.  On 
each  electrode  side  there  is  one  end  plate  with  gas  manifolds 
made  of  Plexiglas15  with  thickness  of  1 .0  cm,  current  col¬ 
lector  plate  made  of  expanded  graphite  with  thickness 
0.03  cm,  two  gas  and  electrolyte  distributor  plates  made 
of  expanded  graphite  each  having  thickness  of  0.03  cm,  a 
gas  diffusion  electrode  composed  of  diffusion  layer  made  of 
carbon  paper  (thickness  0.026  cm)  and  of  catalytic  layer 
made  of  hydrophobized  Pt/C  porous  catalyst  (thickness 
0.01  cm),  and  an  insulating  porous  matrix  with  liquid  elec¬ 
trolyte  hold-up  having  thickness  of  0.026  cm.  Superposition 
of  these  perforated  layers  creates  two  three-dimensional 
systems  of  channels  for  the  flow  of  reactant  gases  and 
electrolyte,  respectively.  The  overall  surface  of  the  PWAFC 
was  196  cm2.  The  active  surface  of  the  electrodes  was 
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Fig.  1.  The  elements  of  the  PWAFC  (physical  model). 
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Fig.  2.  Cross-sectional  view  of  PWAFC.  C i  and  A i  are  designations  for  cathode  and  anode  side  layers;  Taj  and  TAij  are  designations  for  interfacial 
temperatures;  the  thickness  of  layers  is  given  in  mm. 


36  cm2  on  each  side.  The  expanded  view  of  the  PWAFC  is 
shown  in  Fig.  1.  The  cross-sectional  view  with  the  designa¬ 
tion  of  layers,  their  thickness  and  designation  of  interfacial 
temperatures  is  shown  in  Fig.  2.  The  thickness  of  the 
complete  monocell  between  the  end  plates  is  2.66  mm. 
The  disposition  of  channels  for  electrolyte  and  gas  flows 
on  the  distributor  2  at  the  cathodic  side  (see  Fig.  1)  is 
presented  in  Fig.  3  in  order  to  facilitate  the  notion  about 
reactant  gases  and  electrolyte  flow  pattern  in  the  through- 
flow  configuration  of  the  PWAFC.  In  this  configuration  the 
concentrated  PWA  solution  is  recirculated  with  low  flow  rate 
(1.0  cm3  min-1)  through  the  cell  in  direction  from  cathode 
to  anode  side. 

2.2.  PWAFC  structure  and  materials 

The  end  plates  of  the  fuel  cell  were  made  of  Plexiglas.  The 
current  collectors  were  made  of  goldened  copper  plates 
having  thickness  0.25  mm  and  were  placed  into  the  cut  area 
of  the  C5  and  A5  collectors  (see  Fig.  1).  The  distributors 
were  made  of  expanded  graphite  sheets  (Papyex®,  Le 
Carbonne  Loraine,  France).  This  material  is  anisotropic  in 
thermal  and  electric  conductivity:  the  thermal  conductivity 
parallel  to  surface  is  Ax, y  =  11  1  W  m_1  K  1  and  perpen¬ 
dicular  to  the  surface  is  Az  =  0.87  W  m_1  K-1;  the  electric 
conductivity  parallel  to  the  surface  is  <Jx,y  —  1111  S  cm-1 


Fig.  3.  The  disposition  of  channels  for  electrolyte  and  gas  flows  on 
distributor  2  layer. 


and  perpendicular  to  surface  is  ffZ  =  4  S  cm  1  [23]2.  The 
carbon  paper  Toray  TGP  090  having  thickness  of  0.27  mm, 
porosity  cps  =  0.76  and  thermal  conductivity 
As  =  2.508  W  m-1  K  1  was  used  as  a  diffusion  layer  for 
gas  diffusion  electrodes  preparation.  The  electric  conduc¬ 
tivity  of  this  carbon  paper  is  parallel  to  surface 
Ax, y  =  200  S  cm-1  and  perpendicular  to  surface 
Az  1  14.3  S  cm  1  [24].  The  0.1mm  thin  catalyst  layer 
was  screen-printed  on  this  carbon  paper.  The  catalyst  layer 
consisted  of  20  wt.%  Pt  on  active  carbon  hydrophobized 
with  PTFE.  The  final  Pt  load  in  the  electrode  was 
0.5  mg  cm-2  [12].  The  catalyst  layer  is  composed  predo¬ 
minantly  of  porous  active  carbon  for  which  we  take  the 
thermal  and  electric  conductivity  data  of  graphite:  thermal 
conductivity  in  the  range  between  112  and  188  W  nT  1  K  1 
and  electric  conductivity  15,385  S  cm-1  [25], 

2.3.  Description  of  fluid  flow  processes  and  conditions  in 
the  PWAFC 

The  reactant  gases,  hydrogen  (H2)  and  oxygen  (02)  are 
fed  separately  to  the  anode  and  cathode  of  the  PWAFC, 
respectively.  The  gases  first  pass  through  the  tubular  chan¬ 
nels  of  3  mm  i.d.  in  the  end  plates,  then  enter  the  right  lower 
corner  of  each  side  into  the  current  collector  horizontal 
channel  having  a  cross-section  of  3  mm  x  0.3  mm.  They 
pass  then  through  the  horizontal  perforation  of  the  same 
cross-section  in  the  first  distributor  plate  and  enter  into  the 
system  of  parallel  vertical  channels  of  the  second  distributor 
plate.  The  reactant  gases  flow  thus  along  the  diffusion  layer 
of  each  gas  diffusion  electrode,  diffuse  freely  into  the 
catalytic  layer  and  the  gases  in  excess  flow  out  at  the  upper 
left  corner  of  each  electrode  side  in  the  same  manner  as  they 
entered  the  fuel  cell,  but  in  opposite  direction. 

The  flow  of  liquid  electrolyte  (PWA)  is  in  the  direction 
from  cathode  to  anode.  The  electrolyte  enters  the  system  of 
channels,  made  in  essentially  the  same  way  as  those  for 


2  Alternatively,  data  for  pure  expanded  flexible  mineral  graphite  plain 
homogeneous  gasketing  sheets  of  Thermodynamics  Asia  Pacific  Pte  Ltd., 
Singapore  (www.thermody.com/style7000pgl.htm)  or  Le  Carbonne 
Loraine  —  Papyex®  sheet  characteristics. 
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reactant  gases,  at  the  left  lower  comer  of  the  cathode.  When 
the  electrolyte  flow  reaches  the  gas  diffusion  cathode,  it 
flows  through  the  channels  directly  into  the  porous  inert 
matrix.  Here  it  distributes  by  capillary  forces  within  the 
matrix  and  within  the  catalytic  layer  of  both  electrodes.  The 
electrolyte  comes  out  at  the  right  upper  corner  of  the  anode. 
The  flow  of  electrolyte  through  the  porous  inert  matrix  is 
limited.  Therefore,  the  electrolyte  in  excess  is  allowed  to 
come  out  at  the  upper  left  corner  of  the  cathode. 

The  two  systems  of  channels  are  completely  isolated.  The 
reactant  gases  and  electrolyte  come  together  only  in  the 
catalyst  layer  of  both  gas  diffusion  electrodes,  where  they 
form  a  subtly  equilibrated  three-phase  contact  in  which  the 
electrochemical  reactions  proceed. 

The  maximum  work  output  in  the  form  of  electrical 
energy  of  the  ideal  H2/O2  fuel  cell  operating  at  constant 
temperature  and  pressure  is  given  by  the  change  of  the  Gibbs 
free  energy  of  the  overall  reaction  between  H2  and  CK: 

We/dx  =  A  Gr  =  EimxnN0e  (1) 

where  Em ax  is  the  reversible  cell  potential,  n  the  number  of 
moles  of  electrons  transferred  per  mole  of  fuel  consumed 
(n  =  2  in  this  case).  No  =  6.023  x  1023  e  mol  1  (Avogadro’s 
number),  and  e  =  1.6021  x  10  19  C  e-1  is  the  electric  charge 
of  electron.  The  product  Noe  =  F  =  96,487  C  mol-1  e-1 
is  the  Faraday’s  constant.  Then 

£max  =  ~(^jF  =  ]'229Y  (2) 


at  298  K  and  1  bar.  The  reversible  cell  potential  is  given  by 
the  Nernst  equation: 


E  =  E0  + 


f  flH2(flQ2)1/2\ 

y  «h2o  J 


(3) 


It  is  obvious  from  Eq.  (3)  that  the  cell  potential  increases 
with  an  increase  in  the  activity  of  reactant  gases  and  with  a 
decrease  in  the  activity  of  products.  When  the  activity  of 
gases  is  less  than  1,  the  second  term  in  Eq.  (3)  is  negative. 
This  means  that  with  the  increase  of  temperature  the  equili¬ 
brium  cell  potential  will  decrease  and  so  will  the  maximum 
electrical  work  output.  Therefore,  if  the  electrochemical 
reaction  kinetics  permits,  one  should  be  interested  to  work 
at  low  temperatures  in  order  to  obtain  more  electrical  work 
output.  However,  in  practice  the  fuel  cell  voltage  is  much 
lower  than  the  equilibrium  Nernst  potential  since  there  are 
several  irreversible  losses  connected  with  the  existence  of 
the  over-potentials.  The  fuel  cell  voltage  can  thus  be 
expressed  with  the  following  equation: 

^  C  f/act  7/diff  flcond 

where  ?/act  is  the  activation  over-potential,  >7diff  the  diffusion 
over-potential  and  »/con d  is  the  cell  conductivity  over-potential. 
These  sources  of  the  losses  in  the  fuel  cell  voltage  are  not  the 
only  ones.  There  are  several  other  losses  on  the  level  of  the 
fuel  cell  system.  All  these  losses  can  be  expressed  through 


Fig.  4.  Polarization  curve  for  the  PWAFC  at  298  K  and  1  bar  Ft2-02  gas 
pressure  and  the  power  density  as  a  function  of  current  density  [1], 


different  efficiencies  [26].  Three  of  them  are  of  utmost 
importance  for  the  practical  fuel  cell  development: 

•  thermal  efficiency  e*  =  AGr/ AHr  (Gibbs  free  energy 

change/enthalpy  change); 

•  voltage  efficiency  sv  =  V/E  (actual  cell  voltage/thermo¬ 
dynamic  cell  voltage); 

•  current  efficiency  £/  =  efU  (actual  current/current 

equivalent  of  reactant  molar  flow  rate), 

where  eF  is  the  current  efficiency  at  U  =  1,  and  U  is  fuel 
utilization  factor  ( U  =  0  :  1).  The  electrochemical  effi¬ 
ciency  of  the  fuel  cell  is  then  given  as 

eE  =  SthSySF  u  (5) 

For  the  numerical  simulation  of  the  fluid  dynamics  and 
heat  transfer  we  have  used  the  experimental  data  measured 
at  the  peak  power  conditions  of  the  PWAFC  [  1  ]  as  the  input 
values  shown  in  Fig.  4.  The  values  of  these  parameters  are 
current  density  (i  =  1.5  A  cm-2),  cell  voltage  (U  =  0.47  V) 
and  electric  power  density  output  ( Pe\  =  0.7  W  cm-2).  The 
fuel  utilization  in  the  numerical  simulation  was  considered 
to  be  1.0  (U  =  1.0).  It  was  assumed  that  fuel  and  oxygen  are 
consumed  at  their  stoichiometric  ratio  to  form  water.  There¬ 
fore,  £/  =  EpU  =  1.0.  The  theoretical  maximum  reversible 
potential  of  the  H2/O2  fuel  cell  at  room  temperature  (300  K) 
and  normal  pressure  (1  bar)  is  Uq  =  1.229  V  at  the  theore¬ 
tical  maximum  value  of  the  free  enthalpy  change 
A G®  =  —237.2  kJ  mol-1,  and  the  reaction  enthalpy  change 
AH®  =  —285.8  kJ  mol-1,  that  is  equal  to  the  high  heating 
value  (HHV)  of  the  fuel  [26].  The  experimentally  obtained 
voltage  at  peak  power  density  was  0.47  V.  From  these 
values  we  obtain  for  voltage  efficiency  Sv  =  V/E  — 
0.47/1.229  =  0.38,  and  for  the  thermal  efficiency 
£th  =  AGr/A//r  =  (— 237. 2)/(— 285.8)  =  0.83.  The  elec¬ 
trochemical  efficiency  of  PWAFC  at  peak  power  conditions 
is  then  calculated,  according  to  Eq.  (8)  to  be  £g  =  0.317  or 
31.7%. 

Taking  into  account  the  efficiency  at  peak  power  condi¬ 
tions,  we  can  calculate  molar  flow  rates  for  fuel  and 
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oxidant: 


«H 


APd 


A//r°£E 

36.34  x  0.7 
-285.8  x  103  x  0.317 


=  2.80  x  10~4  mol  s_1 


ho  =  =  1-40  x  10  4  mols  1 


(6) 

(7) 


The  corresponding  mass  and  volume  flow  rates  for  hydrogen 
and  oxygen  at  pressure  p  =  1  bar  are  mu  =  5.64  x  10  4 
gs_1,  VH  =  6.99cm3 s_1,mo  =  4.48  x  10~3gs~’,  andl/o  = 
3.49  cm3  s_1,  respectively. 

Before  numerical  modeling,  flow  conditions  in  fuel  cell 
have  to  be  determined.  Since  the  geometry  of  PWAFC  is 
known  in  advance  (see  Section  2.1),  we  have  checked  for 
critical  cross-sections  in  the  system  of  channels  for  gas  flow. 
The  narrowest  part  of  the  hydrogen  and  oxygen  inlet  chan¬ 
nels  has  a  rectangular  cross-section:  width  3  mm  and  height 
0.3  mm.  The  hydraulic  diameter  for  this  cross-section  is 
calculated  /)h  =  0.545  mm.  The  average  linear  velocities  of 
hydrogen  and  oxygen  through  this  cross  section  are 
vh  =  7.77  and  3.88  m  s  ,  respectively.  The  Reynolds  num¬ 
ber  for  the  flow  of  hydrogen  would  then  be 


Re  h 


vh 


7.77  x  0.545  x  10"3 
111  x  10-6 


=  38  «  2320 


(8) 


and  for  oxygen  Reo  =  131  2320.  It  is  obvious  that  even  at 

the  highest  flow  rate  through  the  critical  cross-section  the 
Reynolds  number  is  in  both  cases  deep  in  laminar  flow  type 
region,  therefore  these  findings  are  applied  in  numerical 
model  solution  procedure.  In  this  evaluation  we  have  pre¬ 
supposed  that  the  inlet  gas  is  not  consumed  in  the  electro¬ 
chemical  reaction  (U  =  0),  so  that  the  gas  flow  is  maximal. 
In  practice  the  fuel  utilization  factor  was  U  =  0.66. 

After  describing  the  flow  of  reactant  gases  in  the  PWAFC, 
we  shall  describe  now  the  flow  characteristics  of  the  liquid 
electrolyte  —  a  60  wt.%  water  solution  of  H3PW12O40.  In 
our  through-flow  configuration  of  the  PWAFC  the  concen¬ 
trated  electrolyte  enters  the  fuel  cell  on  the  cathode  side, 
flows  through  the  cell,  as  schematically  shown  in  Fig.  5,  and 
comes  out  diluted  with  the  electrochemically  formed  water 
on  the  anode  side.  The  concentrated  electrolyte  was  pumped 
through  the  cell  with  a  flow  rate  of  1  cm3  min“ 1 .  The  molar 


and  the  mass  flow  rates  of  the  electrochemically-produced 
water  are  calculated  from  the  molar  flow  rates  of  the 
reactants.  The  molar  flow  rate  of  water  is  equal  to  the  molar 
flow  rate  of  hydrogen,  when  the  fuel  utilization  factor 
U=  1.0. 

The  mass  and  volume  flow  rates  of  produced  water  are 

mw  =  Ah  +  mo  =  «hMh2o  =  2.80  x  10~4  x  18 

=  5.04  x  lO^gs'1  (9) 

and  Vw  =  5.05  x  1 0  3  cm3  s  1 .  The  mass  and  volume  flow 
rates  of  the  electrolyte  are  calculated  as  follows: 

VE  =  1  cm3  min~ 1  =  1.66  x  10'2cm3  s_1  (10) 

mE  =  VePe  =  1 .66  x  10~2  x  2.0  =  3.33  x  10~2  g  s_1 

(ID 

where  pE  =  2.0  g  cm-3  is  the  density  of  a  60  wt.%  PWA 
solution  at  298  K  [7],  The  specific  heat  of  the  electrolyte 
solution  is  cE  =  3.4  kJ  kg-1  K  1  and  its  thermal  conductiv¬ 
ity  is  AE  =  0.7  W  m-1  K-1. 

2.4.  Description  of  the  heat  transfer  processes  and 
conditions  in  the  PWAFC 

In  the  heat  transfer  part  of  the  model  we  encounter  several 
different  types  of  heat  transfer  mechanisms  such  as  conduc¬ 
tion,  convection,  diffusion  and  radiation.  We  also  encounter 
several  types  of  material  secondary  structure:  porous  and 
non-porous  solids,  liquids,  electric  conductors  and  insula¬ 
tors.  The  most  complex  is  the  heat  transfer  through  porous 
medium  like  gas  diffusion  electrode.  In  our  case  the  gas 
diffusion  electrodes  were  made  of  0.3  mm  thin  carbon  paper 
with  porosity  <p  =  0.76  on  which  the  0.1  mm  thin  porous 
catalyst  layer  was  screen-printed.  The  heat  conductivity  of 
porous  medium  can  be  estimated  as  the  geometric  mean 
value  of  the  heat  conductivity  of  solid  (As)  and  of  the  liquid 
filling  the  pores  (/.(): 

AG  =  As1-nf  (12) 

This  approximation  holds  when  the  Af  is  not  much  higher 
than  As.  The  heat  conductivity  values  for  hydrogen,  oxygen 
and  for  the  carbon  paper  are  0.183,  0.0268  and 


Fig.  5.  Schematic  representation  of  the  electrolyte  flow  within  PWAFC. 
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0.180  W  m-1  K_1,  respectively.  We  calculate  the  heat  con¬ 
ductivity  for  anode  Aq  =  0.182  Wm  1  K  1  and  the  heat 
conductivity  for  cathode  Ac,  =  0.0423  W  m  1  K-1.  One  can 
notice  that  the  heat  conductivity  of  hydrogen  and  that  of 
anode  material  are  very  similar,  therefore,  the  real  structure 
of  anode  material  is  not  important. 

The  heat  flow  which  is  generated  during  the  fuel  cell 
operation  can  be  calculated  from  the  peak  power  density  and 
the  electrochemical  efficiency  of  the  PWAFC: 

.„  =  Pe i(l  -  fid)  _  0.7  x  10-4  x  (1  -  0.317) 
qg  ~  Eel  “  0.317 

=  1.51  x  10~4WnW2  =  1.51  Went'2  (13) 

This  heat  is  generated  as  a  Joule  heat  due  to  electrolyte 
resistance  and  as  the  electrochemical  reaction  heat  in  the 
three-phase  region  of  the  cathode  catalyst  layer.  To  calculate 
the  specific  volume  heat  generation,  which  one  needs  to 
define  boundary  conditions  for  the  solution  of  numerical 
model,  one  must  first  calculate  the  active  volume  in  which 
the  heat  is  generated.  The  active  surface  of  the  anode  is 
36.34  cm2  and  the  thickness  of  the  inert  matrix  with  liquid 
electrolyte  holdup  is  0.26  mm.  The  thickness  of  the  three- 
phase  electrochemical  reaction  region  is  in  the  order  of 
10-5  mm  and  can  be  neglected.  Thus,  the  active  volume 
for  heat  generation  is  VA  =  0.9448  cm3,  the  specific  volume 
heat  generation  is  q"g  =  55.08  W  crrT3  and  the  absolute 
value  of  the  produced  heat  flow  is  Qg  =  54.87  W. 

The  diffusion  heat  flow  is  generated  when  the  laminar 
flow  of  reactant  gas  (FF  or  O2)  passes  along  the  outer  surface 
of  the  GDE  diffusion  layer  while  the  reactant  gas  is  con¬ 
sumed  in  the  electrochemical  reaction  occurring  within  the 
catalyst  layer  of  the  GDE  which  adheres  to  the  diffusion 
layer. 

In  the  case  of  our  model,  the  fuel  utilization  factor  is 
(7=1  and  the  diffusion  flows  of  the  reactant  gases  in  both 
electrodes  are  equal  to  their  mass  flows.  We  calculated  these 
mass  flows  from  the  peak  power  conditions  of  the  PWAFC. 
From  these  data  we  calculate  the  inlet  enthalpy  flows  of  both 
reactants:  Hu  =  /W|-K'p  i[7H  and  Ho  =  ihocp.oTo- 

The  heat  transfer  occurs  also  by  the  flow  of  liquid 
electrolyte  through  the  fuel  cell.  Although  the  mass  flow 
of  liquid  electrolyte  is  low,  its  specific  heat  is  much 
higher  than  that  for  gases.  The  inlet  enthalpy  flow  of 
electrolyte  is 

/7c,i  =  mcccTc-,  (14) 

and  the  outlet  enthalpy  flow  is 

H e,o  :  7(n  “t“  U/w('w 7\v.o  (15) 

One  of  the  mechanisms  for  heat  transfer  in  the  PWAFC  is 
also  conduction  in  the  expanded  graphite  sheets.  This  mate¬ 
rial  exhibits  anisotropy  in  both  electric  and  heat  conductiv¬ 
ity.  In  such  a  case,  a  three-dimensional  model  for  stationary 
heat  conduction  in  the  anisotropic  medium  is  used  (Fourier 
law).  Since  for  expanded  graphite  Ax  =  Ay  =  11.1  Wm 


K  1  8>  Az  =  0.87  W  nT1  KT1  one  can  use  the  simplified 
three-dimensional  model  in  which  the  temperature  in  lateral 
directions  changes  faster  than  in  perpendicular  direction  of 
the  expanded  graphite  sheet,  in  accordance  with  difference 
in  heat  conductivity. 

The  linear  velocities  of  electrolyte  through  cathode  and 
anode  are  calculated  from  the  electrolyte  flow  rate  through 
the  electrolyte  channel  cross-section  and  they  are  equal  to 
ve,-,c  =  0.267  mm  s-1  and  ve,^A  =  0.352  mm  s_1,  respec¬ 
tively.  These  linear  velocities  of  electrolyte  are  low,  there¬ 
fore,  we  shall  pre-suppose  that  we  have  a  standstill  layer  of 
electrolyte,  having  heat  conductivity  Ae  =  0.7  W  m-1  K  1 
and  we  shall  treat  the  heat  transfer  in  the  electrolyte  as  a  pure 
heat  conduction. 

The  heat  transfer  from  the  surface  of  the  electrode  to  the 
gas  stream  is  calculated  using  the  following  procedure.  The 
laminar  flow  in  the  channel  of  rectangular  cross-section  is 
considered  to  be  fully  developed.  For  the  rectangular  cross- 
section  with  sides  ratio  5:1  (width  of  the  channel 
a  =  1.5  mm,  height  of  the  channel  b  =  0.3  mm)  and  the 
hydraulic  diameter  D h  =  0.5  mm  the  value  of  the  Nusselt 
criteria  would  be  Nuo  =  5.75  [23].  The  heat  transfer  coeffi¬ 
cient  is  then  calculated  for  anode  (H)  and  cathode  (O)  with 
the  following  equation: 


aH 


NudAh 

Dh 


5.75  x  0.183 
5  x  10-4 


2104.5  W  m~2  K_1 


(16) 


and  similarly  ao  =  308.2  W  m~2  K-1. 

The  heat  flow  with  transfer  from  backside  of  electrode  to 
the  gas  stream  is  defined  as 


(\q/)  =  oApC(7c23  -  7o)  for  the  cathode  (17) 


and 


(\q/)  =  aApA(rA23  -  TH)  fortheanode  (18) 

where,  Tc 23  (K)  and  TA2 3  (K)  are  the  temperatures  at  the 
interfaces  schematically  presented  in  Fig.  6. 

The  radiative  heat  transfer  in  the  rectangular 
(1.5  mm  x  0.3  mm)  channels  with  view  factor  0.01  between 
horizontal  and  vertical  sides  is  so  low  that  one  can  neglect  it 
and  take  into  consideration  only  radiative  heat  transfer 
between  two  parallel  surfaces  of  layers: 

k  =  «■*■(,/„)  +  (I/t2>  -.(fl.-AO  <»> 


Emitivities  for  electrode  (carbon  paper,  £1  =  0.95)  and 
distributor  (expanded  graphite,  £2  =  0.95)  materials,  the 
value  of  Stefan-Boltzmann  constant  (cR  =  5.67  x  1 0  6 
W  m-2  K  4)  and  r34  (K)  is  either  7b34  or  7’A34  (see 
Fig.  6),  and  723  (K)  is  either  7c23  or  rA23  (see  Fig.  6). 

Although  we  have  developed  the  mathematical  model  of 
heat  transfer  in  the  PWAFC  we  decided  to  find  numerical 
rather  than  analytical  solution,  due  to  the  complex  three- 
dimensional  geometry  of  the  fuel  cell.  For  this  purpose  we 
have  used  the  CFD  software  package  PHOENICS  [22]. 
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TC67  TC56  TC34  TC12  TA12  TA34  TA56  TA67 

TC45  TC23  Tl  TA23  TA45 


Fig.  6.  The  scheme  of  heat  transfer  mechanisms  in  the  PWAFC. 


3.  Numerical  model  of  the  PWAFC 

In  the  previous  section  we  have  described  the  physical 
model  of  the  PWAFC  and  the  processes  and  conditions 
during  its  work  which  enable  the  definition  of  boundary 
conditions  for  numerical  modeling  of  mass  and  heat  transfer 
processes  within  the  fuel  cell.  The  scope  of  this  work  was  not 
to  calculate  accurate  numerical  values  of  physical  para¬ 
meters  in  one  point,  but  rather  to  establish  global  relations 
and  distributions  of  particular  parameters  through  the  entire 
fuel  cell. 

There  are  three  main  causes  for  the  heat  generation  within 
the  fuel  cell: 

1.  Electrochemical  (electrocatalyzed)  reaction  between 
mobile  protons  and  adsorbed  oxygen  anions  on  the  Pt 
catalyst  of  cathode. 

2.  Chemical  (heterogeneously  catalyzed)  reaction  between 
molecular  hydrogen,  crossing  over  from  anode  to 
cathode  due  to  too  high  permeability  of  the  membrane 
electrolyte  or  matrix  filled  with  liquid  electrolyte,  and 
adsorbed  oxygen  on  Pt  catalyst  in  the  cathode. 

3.  Ohmic  resistance  of  the  electrolyte  (solid  or  liquid).  The 
resistance  of  ionic  conductors  is  roughly  1000  times 
higher  than  the  resistance  of  electronic  conductors.  The 
heat  generated  in  this  way  is  proportional  to  the  square 
of  current  passing  the  electrolyte  (Joule  heat). 

The  generation  of  heat  within  a  fuel  cell  is  a  direct 
mapping  of  the  distribution  and  consumption  of  reactant 
gases  within  the  porous  gas  diffusion  electrodes.  But,  since 
the  distribution  of  reactant  gases  depends  on  the  fluid  (gases 


and  electrolyte)  dynamics  within  given  fuel  cell  geometry, 
the  study  of  heat  generation  in  the  fuel  cell  is  closely  related 
to  the  study  of  distribution  and  flow  of  fluids  in  it  [19].  The 
mass  flow  of  oxygen  in  the  cathode  and  hydrogen  in  the 
anode  is  directly  proportional  to  the  current  density.  The 
square  of  current  multiplied  by  internal  resistance  of  the  fuel 
cell  (which  is  equal  to  the  resistance  of  ionic  conductor,  i.e. 
electrolyte)  represents  the  Joule  heat. 

Based  on  these  facts  we  have  approached  to  the  study  of 
heat  flow  by  using  the  software  package  for  modeling  the 
fluid  dynamics  and  heat  transfer.  The  results  obtained  with 
numerical  model  on  the  reactant  gas  distribution  in  the  gas 
diffusion  electrodes  (cathode  and  anode)  are  fully  correlated 
with  the  heat  generation  in  the  fuel  cell. 

We  have  also  performed  the  numerical  simulation  of  the 
liquid  electrolyte  flow  through  the  fuel  cell  since  it  regulates 
the  formation  of  the  three-phase  region  within  the  GDE 
catalyst  layer  where  the  electrochemical  reactions  proceed 
and  since  it  regulates  the  liquid  electrolyte  specific  con¬ 
ductance  and  Joule  heat  generation. 

The  combination  of  the  gas  flows  and  electrolyte  flow 
numerical  simulation  within  the  PWAFC  should  give  us  the 
answer  on  the  following  question:  under  which  simulated 
operating  conditions  of  the  PWAFC  one  can  expect  anoma¬ 
lies  that  we  have  encountered  during  experimental  runs  of 
PWAFC. 

3.1.  General  description  of  numerical  modeling  procedure 

The  governing  partial  differential  equations  (PDEs)  are 
solved  for  the  conservation  of  mass,  momentum,  energy  and 
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chemical  species  as  applicable.  The  governing  PDEs,  each 
representing  the  continuity  equation  when  <f>  =  1 ,  is  written 
in  the  form  of  a  Cartesian  tensor  as 

+  V(rp<^v)  -V(r*V(0))=  S  (20) 

yt  convection  term  diffusion  term  source  term 

transfer  term 

where  r  is  the  volume  part,  p  the  density,  </>  the  conserved 
mass,  momentum,  energy  or  species  quantity,  t  the  time,  v 
the  velocity  vector,  %  the  transport  coefficient  (viscosity,  heat 
conductivity,  etc.),  and  S  is  the  change  of  property. 

Each  PDE  is  reduced  within  PHOENICS  software  pack¬ 
age  to  its  finite  volume  form  by  integrating  it  over  the  control 
volumes  of  the  design  domain. 

In  the  numerical  calculation  other  additional  equations  are 
used  which  describe  the  thermodynamic  properties  (density, 
etc.),  transport  properties  (viscosity,  diffusion,  convection, 
etc.),  generating  properties  (radiation  absorption,  chemical 
generation,  viscous  dissipation,  etc.)  and  also  interface 
transport  properties  (moment,  energy,  etc.).  The  boundary 
conditions  have  to  be  set  as  a  constant  value  of  certain 
property,  as  a  constant  flux  of  certain  property  and  as  linearly 
and  non-linearly  defined  properties. 

3.2.  Entiy  data 

We  have  developed  a  three-dimensional  model  of  the  fuel 
cell  in  which  the  heat  transfer  proceeds  by  conduction  and 
convection.  The  radiative  heat  transfer  is  negligibly  small 
and  omitted  from  the  model  because  of  the  low  fuel  cell 
working  temperatures  (less  than  353  K)  and  because  of 
small  temperature  difference  between  the  fuel  cell  layers  in 
the  perpendicular  direction.  Based  on  the  preliminary  estima¬ 
tion  of  the  mass  transfer  conditions  in  the  fuel  cell  (see 
Section  2.3),  we  have  chosen  the  way  of  solving  the  numeri¬ 
cal  model,  which  is  adapted  for  the  laminar  flow  of  the  gases. 

The  whole  model  of  the  fuel  cell  is  oriented  as  shown  in 
Fig.  1,  where  the  in-plane  directions  are  X-  and  F-,  and  the  Z- 
direction  is  perpendicular  to  the  plane.  The  geometrical 
model  for  numerical  simulation  is  composed  of  nine  indi¬ 
vidual  layers,  each  having  the  thickness  of  0.30  mm,  with 
the  exception  of  the  middle  layer  representing  the  inert 
porous  matrix  with  thickness  of  0.26  mm.  The  outer  dimen¬ 
sions  of  the  model  are,  therefore,  140  mm  in  X-direction, 
140  mm  in  T-direction  and  2.66  mm  in  Z-direction. 

The  porous  elements  of  the  fuel  cell,  such  as  both  gas 
diffusion  electrodes  are  set  to  a  value  of  0.76  of  volume 
porosity.  The  numerical  calculation  of  the  model  was 
divided  into  two  parts: 

1.  the  calculation  of  the  oxygen  flow  through  the  gas 

channels  on  the  cathode  side; 

2.  the  calculation  of  the  electrolyte  flow  through  the 

channels  across  the  entire  fuel  cell. 

In  calculations  of  the  first  part  we  defined  the  electrolyte 
channels  on  the  cathode  side  as  filled  blocks  (standstill 


electrolyte).  The  oxygen  flow  for  the  given  operation  con¬ 
ditions  of  the  fuel  cell  is  known.  The  linear  velocities  and 
mass  flow  rates  were  calculated  with  regard  to  the  inlet  and 
outlet  cross-sections,  which  are  24  mm2  (3  mm  x  8  mm). 

The  outlet  cross-section  had  to  be  defined  at  the  0  kPa 
pressure.  The  proper  sources  for  the  heat  generation  were 
defined  and  the  values  of  specific  heat  flow  per  unit  volume 
were  calculated. 

3.3.  Numerical  calculations 

The  structured  grid  was  developed  based  on  the  geome¬ 
trical  model  of  the  fuel  cell.  The  grid  is  composed  of 
130  x  58  x  21  (158340)  nodes.  These  nodes  are  distributed 
across  the  X-F-plane  in  such  a  way  as  to  enable  the  con¬ 
tinuous  representation  of  computed  results  in  all  regions.  In 
the  regions  containing  the  gas  and  electrolyte  distribution 
channels  the  grid  is  denser.  To  reach  the  adequate  accuracy 
limits  in  the  next  calculation  step  about  1000  iterations  were 
necessary,  which  took  more  than  15  hours  of  processing  time 
on  a  PC  with  Pentium  233  MMX. 


4.  Analysis  of  the  numerical  simulation  results 


4.1.  Analysis  of  the  numerical  simulation  results  for  the 
oxygen  flow  in  the  cathode  side  half-cell  of  PW AFC 


The  oxygen  flow  in  the  cathode  side  half-cell  was  calcu¬ 
lated  without  heat  generation.  The  pressure  distribution  and 
the  velocity  vectors  were  determined  without  the  influence 
of  the  temperature.  Four  plains  represent  each  layer  within 
the  cathode  side  half-cell.  The  velocity  vectors  and  their 
components  are  presented  separately  for  the  distributors  (left 
side)  and  for  the  porous  cathode  (right  side)  in  Figs.  7-10. 
The  calculated  values  of  variables  have  the  following  de¬ 
signation: 


PI 

U1 

VI 

W1 

Vector 

TEM1 


relative  pressure  (the  value  of  surrounding 
pressure  is  100  kPa)  (Pa); 
velocity  component  in  the  X-axis  direction 
(m  s"1); 

velocity  component  in  the  F-axis  direction 
(m  s_1); 

velocity  component  in  the  Z-axis  direction 
(ms"1); 

velocity  vector  (m  s 
absolute  temperature  (K). 


The  intensity  for  each  represented  variable  is  given  in  gray 
scale.  In  order  to  improve  graphical  visualization  of  each 
variable  value  the  view  in  Z-axis  direction  is  expanded  50 
times.  The  left  side  of  Figs.  7-10  represents  the  oxygen 
velocity  distribution  within  the  channels  of  the  collector, 
distributors  1  and  2,  and  the  right  side  of  Figs.  7-10 
represents  the  oxygen  velocity  distribution  on  the  backside 
plane  of  the  porous  cathode. 
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Fig.  7.  X-axis  component  of  oxygen  flow  velocity  in  the  channels  formed  by  superposition  of  collector,  distributors  1  and  2  (left)  and  on  the  back  surface  of 
the  cathode  (right). 


Fig.  7  (left  side)  reveals  that  the  X-axis  component  of 
velocity  is  uniform  along  the  height  of  the  channel  (3  mm), 
while  it  changes  drastically  along  the  width  of  the  channel 
(0.6  mm).  On  the  front  wall  of  the  channel  it  reaches  3.3  m 
s~ 1  and  then  it  drops  to  nearly  zero  value  on  the  back  wall  of 
the  channel.  One  can  notice  that  the  velocity  of  oxygen  flow 
along  the  channel  rises  and  drops  intermittently,  according 
to  the  widening  and  narrowing  of  the  channel  formed  by 


superimposed  layers  of  collector  and  distributor  1  (see 
Fig.  1).  The  distribution  of  the  X-axis  velocity  component 
in  the  upper  outlet  channel  (the  oxygen  exits  at  the  upper  left 
corner)  is  nearly  the  same  as  that  in  the  lower  inlet  channel, 
but  rotated  for  1 80°  clockwise  around  the  Z-symmetry  axis. 

The  T-axis  component  of  the  oxygen  flow  velocity  dis¬ 
tribution  is  shown  in  Fig.  8.  In  the  distributor  2  vertical 
channels  the  highest  positive  values  of  this  velocity  compo- 


PH0T0N 


U1 

_-0 . 73 

5-0.59 

5-0 . 46 

“-0.3Z 

L-0.19 

-0.05 

0.08 

0.ZZ 

0.35 

0.49 

0.6Z 

0.76 

0.90 

1.03 

1.17 


2  X 


Fig.  8.  F-axis  component  of  oxygen  flow  velocity  in  the  channels  formed  by  superposition  of  collector,  distributors  1  and  2  (left)  and  on  the  back  surface  of 
the  cathode  (right). 
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Fig.  9.  Z-axis  component  of  oxygen  flow  velocity  in  the  channels  formed  by  superposition  of  collector,  distributors  1  and  2  (left)  and  on  the  back  surface  of 
the  cathode  (right). 


nent,  about  1  ms-1,  are  achieved  at  the  entrance  into  the  first 
vertical  channel  and  at  the  exit  from  the  last  vertical  channel 
(Fig.  8,  left).  The  Taxis  velocity  component  values  on  the 
backside  surface  of  the  cathode  (Fig.  8,  right)  are  the  highest 
and  positive  in  vertical  stripes  which  match  the  position  and 
geometry  of  the  vertical  distribution  channels  in  the  dis¬ 
tributor  2. 


Fig.  9  represents  the  distribution  of  the  Z-axis  gas  velocity 
component  in  the  cathode  compartment  of  the  PWAFC.  One 
can  notice  that  the  non-zero  values  are  located  only  at  the 
points  of  oxygen  inlet  and  outlet  and  the  points  of  the 
vertical  channel  entrance  and  exit  (Fig.  9,  left). 

The  integral  view  of  the  oxygen  flow  in  the  cathode 
compartment  is  presented  in  Fig.  10  in  the  form  of  velocity 


Fig.  10.  Velocity  vectors  of  oxygen  flow  in  the  cathode  compartment  of  PWAFC. 
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Fig.  1 1 .  The  oxygen  pressure  field  within  the  distribution  channels  of  the  cathode  compartment  (left)  and  on 


the  backside  surface  of  the  cathode  (right). 


vectors.  The  maximum  values  of  the  velocity  vectors  within 
the  channels  (Fig.  10,  left)  are  equal  to  the  maximum  values 
of  the  X-axis  velocity  component  (compare  with  Fig.  7,  left). 
The  velocity  vectors  distribution  on  the  backside  surface  of 
the  cathode  (Fig.  10,  right)  matches  well  the  distribution  of 
the  T-axis  velocity  component  (compare  with  Fig.  9,  right). 

The  pressure  distributions  along  the  channels  (Fig.  11, 
left)  and  on  the  backside  cathode  surface  (Fig.  11,  right)  give 
us  the  best  visual  presentation  of  the  stationary  oxygen  flow 
conditions  in  the  cathode  compartment  of  the  PWAFC.  The 
overall  pressure  drop  between  the  oxygen  inlet  point  and  the 
oxygen  outlet  point  amounts  64  Pa  (Fig.  11,  left).  One  can 
notice  that  the  pressure  uniformly  drops  along  the  horizontal 
double-width  channel  from  the  inlet  point  to  the  point  of 
entrance  into  the  first  two  vertical  channels.  Here,  the 
pressure  somewhat  rises  due  to  the  horizontal  channel 
narrowing,  and  then  continues  to  drop  uniformly  until  the 
next  branching  off  into  the  vertical  channels.  It  is  important 
to  note  that  the  pressure  drop  along  the  horizontal  channel  is 
about  twice  as  large  as  the  pressure  drop  along  the  vertical 
channels.  The  stationary  picture  of  pressure  drop  on  the 
backside  surface  of  the  cathode  reflects  the  one  in  the  last 
layer  of  the  distributor  2,  with  the  exception  that  the  pressure 
field  is  spread  over  the  entire  surface  of  the  cathode  (Fig.  1 1, 
right).  One  can  observe  the  stepwise  pressure  drops  between 
the  regions  delineated  by  the  vertical  electrolyte-filled  chan¬ 
nels.  These  channels  are  impermeable  for  the  gas  in  our 
numerical  model.  The  values  of  these  pressure  drops  along 
both  sides  of  the  vertical  electrolyte  filled  channels  are  equal 
to  the  pressure  drops  in  both  horizontal  distribution  channels 
(lower  inlet  and  upper  outlet  channel). 

The  analysis  of  the  pressure  field  reveals  that  the  highest 
pressure  drops  are  generated  within  inlet  and  outlet  hori¬ 


zontal  gas  distribution  channels  formed  by  superposition  of 
the  collector  and  distributor  1.  The  immediate  practical 
conclusion  can  be  drawn:  the  cross-section  of  these  two 
channels  should  be  increased  and  the  obstacles  created  by 
intermittent  narrowing  and  widening  along  these  channels 
should  be  removed.  Also,  the  gas  inlet  and  outlet  should  be 
enabled  on  both  sides  of  the  horizontal  channels. 

4.2.  Analysis  of  the  numerical  simulation  results  for  the 
electrolyte  flow  within  the  PWAFC 

The  electrolyte  flow  within  the  entire  fuel  cell  geometry 
was  numerically  modeled,  but  in  this  case  the  gas  flow  was 
excluded.  The  flow  scheme  of  electrolyte  is  slightly  different 
from  that  of  the  reactant  gases.  This  flow  pattern  is  shown  in 
Fig.  5. 

The  volume  flow  rate  of  the  electrolyte  is  low  and  in  the 
model  we  took  the  value  of  1  cm3  min"  1 .  No  doubt,  this  flow 
is  laminar,  because  the  inlet  linear  velocity  through  the 
9  mm2  cross-sectional  area  is  1.856  x  10~3  m  s-1.  The  X- 
axis  component  of  electrolyte  flow  velocity  distribution 
within  PWAFC  is  shown  in  Fig.  12  (left).  These  velocity 
component  values  are  the  highest  in  the  horizontal  inlet  and 
outlet  channels,  and  are  equal  to  zero  in  the  vertical  channels 
and  close  to  zero  in  the  porous  matrix.  The  highest  values  of 
the  X-axis  velocity  component  do  not  exceed 
2  x  1 0  2  ms-1.  The  Z-axis  component  of  the  electrolyte 
flow  velocity  distribution  is  also  close  to  zero  in  the  channels 
and  in  the  porous  matrix,  as  shown  in  Fig.  12  (right). 
However,  the  most  interesting  is  the  distribution  of  the  Y- 
axis  electrolyte  velocity  component  in  the  vertical  channels 
of  the  distributor  2  and  of  electrodes,  and  in  the  inert  porous 
matrix  with  electrolyte  holdup  (Fig.  13).  In  the  given  upflow 
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Fig.  12.  The  X-axis  velocity  component  (left)  and  the  Z-axis  velocity  component  (right)  of  the  electrolyte  flow  distributions. 


configuration  of  the  electrolyte  circulation  through  the  fuel 
cell,  the  electrolyte  flow  is  divided:  one  part  of  electrolyte 
flows  through  porous  matrix  and  exits  through  the  upper 
horizontal  channel  on  the  anode  side;  the  other  part  of  the 
electrolyte  flows  upward  the  vertical  channels  in  the  cathode 
and  upward  the  porous  matrix  and  then  exits  through  the 
upper  horizontal  channel  on  the  same,  cathode  side. 

As  in  the  case  of  oxygen  flow,  the  most  illustrative 
representation  of  the  stationary  electrolyte  flow  conditions 


within  the  PWAFC  is  obtained  by  examining  the  pressure 
field  (Fig.  14).  The  highest  pressure  is  created  at  the  inlet 
point  of  electrolyte  flow  which  amounts  to  5.4  Pa.  Especially 
informative  is  the  pressure  field  created  in  the  porous  inert 
matrix  which  determines  the  rate  of  electrolyte  renewal 
within  the  three-phase  region  of  the  electrode  catalyst  layer 
that  adheres  to  the  porous  matrix.  This  rate  in  turn  deter¬ 
mines  the  rate  of  the  electrochemical  reactions  on  both 
electrodes. 
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Fig.  13.  The  7-axis  velocity  component  of  the  electrolyte  flow  distribution. 


As  can  be  seen  in  Fig.  14  (right)  the  created  electrolyte 
pressure  field  within  the  porous  matrix  is  far  from  being 
homogeneous.  The  ratio  between  the  pressure  drop  along  the 
horizontal  distribution  channel  and  along  the  vertical  chan¬ 
nels  and  in  vertical  direction  of  the  matrix  is  lower  than  it 
was  in  the  case  of  gas  pressure  drop  in  horizontal  channels 
and  in  the  vertical  direction  of  the  backside  surface  of  the 
cathode. 

When  looking  from  the  cathode  side,  the  upper  right 
corner  of  the  matrix,  and,  on  the  same  place  also  both 
electrodes,  will  be  poorly  wetted,  while  the  diagonally 
positioned  lower  left  corner  of  the  matrix  and  of  both 
electrodes  will  be  wetted  sufficiently.  When  one  makes 
the  superposition  of  the  schematically  represented  pressure 
gradients  of  the  electrolyte  in  matrix  and  of  the  gas  (oxygen) 
in  the  cathode,  the  picture  like  the  one  represented  in  Fig.  15 
emerges. 

It  is  evident  that  the  lower  right  corner  of  the  catalyst  layer 
on  cathode  will  be  well  fed  with  oxygen,  but  will  not  be 


Cathode 


Fig.  15.  Schematic  representation  of  the  pressure  gradients  of  gas  in  the 
electrode  and  of  electrolyte  in  the  matrix. 


wetted  enough  with  fresh  electrolyte.  The  lack  of  electrolyte 
(or  inhomogeneous  electrolyte  distribution)  can  cause  sev¬ 
eral  deteriorating  effects.  In  such  cases  the  electrolyte  would 
be  diluted  faster  with  the  electrochemically  produced  water 
and  the  local  Ohmic  resistance  would  increase.  Conse¬ 
quently,  the  local  temperature  would  also  increase.  This 
would  imply  local  drying  of  the  matrix.  On  the  other  side, 
the  pressure  of  both  gases  is  highest  exactly  at  this  point,  and 
this  creates  favorable  conditions  for  the  hydrogen  crossover 
with  all  possible  deterioration  consequences. 

4.3.  Analysis  of  the  temperature  distribution  in  the  PWAFC 

The  temperature  distribution  in  the  electrolyte  and  in  the 
middle  of  the  matrix  filled  with  electrolyte  was  calculated 
for  the  following  conditions: 

•  the  cell  works  at  the  peak  power  conditions; 

•  the  surrounding  of  the  cell  is  adiabatic; 

•  the  heat  is  generated  uniformly  at  the  cathode-side  surface 
of  the  matrix; 

•  the  cell  is  cooled  only  with  the  flowing-in  electrolyte 
having  the  temperature  of  293  K. 

The  temperature  fields  across  the  cell  and  in  the  central 
layer  of  the  matrix  resemble  somewhat  the  analogous  elec¬ 
trolyte  pressure  fields,  as  one  would  expect  (Fig.  16).  The 
contours  of  the  position  of  the  horizontal  and  vertical 
distribution  channels  for  gas  and  electrolyte  can  be  recog¬ 
nized  in  the  temperature  field  of  the  matrix  (Fig.  16,  right). 
The  temperatures  along  the  vertical  distribution  channels  for 
the  electrolyte  are  somewhat  lower  than  in  the  other  places  in 
the  matrix.  But,  it  is  obvious  that  the  circulation  of  the 
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Fig.  16.  The  temperature  distribution  across  the  cell  (left)  and  in  the  central  matrix  layer  (right). 


electrolyte  alone  does  not  suffice  for  effective  cooling  of  the 
cell  at  the  peak  power  conditions.  In  order  to  improve  the 
cooling  of  the  fuel  cells  within  the  fuel  cell  stack  it  would  be 
necessary  to  introduce  the  heat  exchangers  with  water  as 
cooling  medium.  Only  by  introducing  the  effective  cooling 
system  one  can  assure  the  normal  cell  operating  conditions 
that  would  prevent  the  temperature  to  rise  above  the  water 
boiling  point  at  the  surrounding  pressure. 

5.  Validation  of  the  numerical  model  results  by 
designed  gas  flow  experiment 

The  original  idea  of  numerical  model  validation  was  to 
design  experiment  in  which  the  basic  parameters  of  the  fluid 
flow  (the  mass  flows  of  reactant  gases  and  of  the  electrolyte) 
and  of  the  temperature  distribution  within  the  fuel  cell  would 
be  measured  continuously  in  PWAFC  with  given  geometry 
and  operational  conditions.  This  would  demand  the  mea¬ 
surement  of  temperature  in  the  catalyst  layer  of  the  cathode 
and  on  the  surface  of  other  layers  inside  the  PWAFC  on  at 
least  20  points.  Such  an  experiment  design  brings  about 
some  inherent  drawbacks  like  substantially  diminishing  the 
active  surface  of  the  cathode  catalyst  layer  due  to  tempera¬ 
ture  sensors  positioning,  corrosion  of  the  sensors  and  wiring, 
problems  in  assuring  the  tightness  between  layers  containing 
temperature  sensors,  etc. 

We  decided,  therefore,  to  simulate  the  development  of  the 
stationary  flow  conditions  of  oxygen  in  the  cathode  com¬ 
partment  of  the  PWAFC  by  assembling  the  cathode-side 
half-cell,  putting  the  sheet  of  the  whiteprint  paper  on  the 
catalyst  layer,  covering  it  with  a  transparent  Plexiglas  “ 
plate,  and  recording  the  time  evolution  of  the  coloring 
reaction  of  diazo  dye  when  exposed  to  white  light  and  a 


flow  of  ammonia.  Since  the  flow  characteristics  of  ammonia 
are  close  to  those  of  oxygen  and  the  whiteprint  paper 
coloring  reaction  is  fast,  this  experiment  offers  the  fluid 
flow  and  chemical  reaction  conditions  that  approach  most 
closely  to  the  conditions  of  oxygen  flow  and  the  electro¬ 
reduction  of  oxygen  in  the  cathode  compartment.  On  the 
other  side  the  lateral  temperature  distribution  in  the  cathode 
depends  on  the  Joule  heat  generated  within  the  fuel  cell.  The 
latter  is  proportional  to  the  square  of  current  density  lateral 
distribution,  which  in  its  turn  is  directly  proportional  to  the 
lateral  distribution  of  oxygen  mass  flow  within  the  porous 
medium  of  gas  diffusion  cathode.  The  dimensions  of  the 
particular  layers  in  this  half-cell  are  the  same  as  those  used  in 
the  geometry  entry  of  the  numerical  model. 

The  whiteprint  paper  was  mounted  on  the  catalyst  layer 
surface  of  the  cathode  and  the  half-cell  was  then  closed  by 
another  transparent  Plexiglas®  end  plate.  The  time  evolution 
of  the  colored  zone  was  filmed  with  the  digital  camera.  The 
whole  experimental  installation  is  schematically  represented 
in  Fig.  17. 

Simultaneously  with  the  video  recording  and  scanning  of 
the  time  evolution  of  gas  distribution  on  cathode,  the  flow 
rate  of  ammonia  was  measured  with  the  soap  bubble  meter 
flow  calibrator  and  the  pressure  drop  of  ammonia  in  the 
cathode  compartment  was  measured  with  the  U-manometer 
filled  with  Squalan  (p  =  0.808  g  cm-3).  The  experiments 
were  performed  at  three  values  of  ammonia  flow  rate:  5.5, 
24.8,  and  49.4  cm3  min-1  which  are  close  to  the  nominal 
values  used  in  the  ammonia  flow  numerical  simulations  (5, 
25  and  50  cm3  min  ~ 1 ).  The  estimated  ammonia  flow  rate 
measurement  uncertainty  is  around  1.1  cm3  min-1.  The 
experimentally  determined  static  pressure  drops  at  three 
different  ammonia  flow  rates  were  8.9,  38.0  and  66.3  Pa. 
The  estimated  pressure  drop  measurement  uncertainty  was 


I.  Lavric  et  al.  /  Journal  of  Power  Sources  96  (2001)  303-320 


317 


nh3 


n2 


Fig.  17.  The  scheme  of  the  experimental  setup  for  measurements  of  the  whiteprint  paper  coloration  with  ammonia. 


4.04  Pa.  In  order  to  compare  the  experimental  results  with 
the  results  of  numerical  simulation,  we  have  calculated  also 
the  dynamic  pressure  for  the  highest  ammonia  flow  rate 
Pd  =  1.5  x  10  3  Pa.  The  dynamic  pressure  at  the  inlet  and 
at  the  outlet  of  the  cathode  compartment  is  low  and  we  shall 
neglect  it  in  our  further  discussion. 

The  video  sequences  acquisition  was  performed  with 
FPS60  video  card  (Fast  Multimedia,  Germany)  and  Premiere 
software  package  (Adobe)  [27].  The  frequency  of  acquisi¬ 
tion  was  one  snapshot  per  second. 

The  pictures  were  then  processed  further  on  with  a  soft¬ 
ware  package  Picture  Publisher  6.0  (Micrografx,  USA)  [28]. 
The  contours  obtained  by  processing  the  successive  snap¬ 
shots  were  then  superposed  in  order  to  obtain  one  picture 
demonstrating  the  time  evolution  of  the  reaction  front. 

The  measured  and  calculated  flow  characteristics  of 
ammonia  are  linear  (Fig.  18).  Using  the  Darcy’s  equation 
for  pressure  drop  for  laminar  flow  we  get  the  following 
equation  of  flow  characteristic  in  the  cathode  compartment 
of  the  fuel  cell,  giving  the  linear  relation  between  the  gas 
velocity  and  pressure  drop: 


Fig.  18.  Measured  and  calculated  ammonia  flow  characteristics. 


Ap  =  KvpV  =  98, 676  x  vpV  (Pa)  (21) 

where  for  ammonia  v  =  14.7  x  1CD6  m2  s-1  and 
p  =  0.6894  kg  m  \  One  can  notice  that  the  measured 
and  calculated  values  for  pressure  drop  agree  well  in  case 
I,  that  is  at  low  ammonia  flow  rate.  At  higher  flow  rates  the 
measured  and  calculated  flow  characteristics  part  and  the 
calculated  values  are  higher  than  measured  ones.  We  ascribe 
this  to  enhanced  leakage  of  ammonia  from  experimental 
assembly  at  higher  pressure  drops. 

The  time  evolution  of  the  reaction  front  measured  at  three 
different  flow  rates  of  ammonia  through  cathode  compart¬ 
ment  of  the  PWAFC  is  represented  in  Fig.  19.  The  coloring 
reaction  front  propagates  faster  along  the  first  vertical 
channel  (on  the  right  side)  and  with  time  lags  needed  for 
the  propagation  along  the  horizontal  channel  (lower  one)  it 
then  propagates  also  along  the  successive  vertical  channels 
(Fig.  19,  case  I).  At  higher  ammonia  flow  rates  the  propaga¬ 
tion  of  coloring  reaction  along  the  first  vertical  channel 
becomes  increasingly  faster  comparing  to  the  propagation 
along  successive  vertical  channels.  Ammonia  from  the  first 
vertical  channel  reaches  the  upper  horizontal  channel  and 
spreads  along  it  faster  than  from  successive  vertical  chan¬ 
nels.  Moreover,  one  can  notice  that  the  ammonia  (and 
consequently  the  coloring  reaction)  propagates  from  the 
upper  horizontal  channel  into  the  last  two  vertical  channels 
downward!  This  creates  the  asymmetry  in  the  propagation 
front  of  the  coloring  reaction  visible  in  Fig.  19  (cases  II  and 
III).  The  propagation  in  the  first  vertical  channel  of  case  II  is 
equal  to  the  propagation  in  the  second  vertical  channel  of 
case  III,  since  the  flow  rate  in  case  III  is  doubled  with  respect 
to  flow  rate  in  case  II.  Detailed  inspection  reveals  that  the 
difference  in  the  reaction  propagation  velocity  along  the  first 
and  second  vertical  channels  increases  with  flow  rate  when 
the  flow  rates  are  high  enough.  This  effect  is  represented  in 
Fig.  20. 
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Fig.  19.  The  time  evolution  of  the  coloring  reaction  front  on  whiteprint  paper  at  three  values  of  ammonia  flow  rate  through  the  cathode  compartment: 
5.5  cm3  muT1  (I);  24.8  cm3  min-1  (II);  49.4  cm3  min-1  (III). 
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Fig.  20.  The  difference  in  reaction  propagation  between  first  and  second 
vertical  channels  for  different  flow  rates  of  ammonia. 


Figs.  19  and  20  demonstrate  the  break  of  diagonal  sym¬ 
metry  in  propagation  of  the  coloring  reaction  zone  when  the 
ammonia  flow  rate  is  increased  from  5.5  to  49.4  cm3  min~ 1 . 
With  the  increase  of  the  ammonia  flow  rate  the  propagation 
of  the  reaction  zone  along  the  first  vertical  channel  is 
becoming  much  faster  than  the  propagation  of  the  reaction 
zone  along  successive  vertical  channels. 


6.  Analysis  of  heat  generation  within  PWAFC 

At  low  reactant  gas  flow  rates  the  gas  is  uniformly 
distributed  over  the  entire  electrode  surface.  The  current 
density  is  low  and  so  is  the  Joule  heat  generated  by  electro¬ 
lyte  resistivity.  The  flow  rate  of  the  electrolyte  is  high 
enough  to  wet  the  entire  inert  matrix.  The  pressure  drop 


Fig.  21.  A  critical  region  of  high  heat  generation  in  the  PWAFC  (ellipse)  on  cathode  (left)  and  on  matrix  (right). 
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Fig.  22.  Superposition  of  the  photocopied  damaged  zone  after  post-mortem  PWAFC  analysis  on  the  gas  pressure  distribution  on  the  backside  of  the  cathode 
obtained  with  the  numerical  simulation  of  the  PWAFC. 


due  to  reactant  gas  flow  is  lower  than  the  bubble  pressure  of 
the  electrolyte  in  the  matrix  and  no  reactant  gas  crossover 
occurs. 

When  the  reactant  gas  flow  rate  increases,  the  heat  load  in 
the  fuel  cell  increases  due  to  higher  current  densities 
involved  and  due  to  higher  rate  of  Joule  heat  generated. 
Moreover,  the  gas  pressure  drops  become  higher  and  the  gas 
is  non-uniformly  distributed  over  the  electrode  surface.  If 
the  local  mismatch  between  the  gas  and  electrolyte  flows  is 
created,  the  locally  generated  heat  can  cause  the  electrolyte 
to  dry.  The  solid  crystalline  electrolyte  is  permeable  for  the 
reactant  gas  (especially  hydrogen)  and  crossover  becomes 
possible.  The  hot  spots  are  created  within  the  cathode 
catalyst  layer  that  can  locally  melt  the  PTFE  and  burn  the 
catalyst.  This  would  lead  to  the  lower  fuel  cell  efficiency  and 
finally  to  the  fuel  cell  break  down.  Based  on  the  results  of  the 
numerical  simulation  of  the  reactant  gas  and  electrolyte 
flows  explained  in  Sections  4.1  and  4.2  and  according  to 
Fig.  15,  one  can  expect  that  the  critical  region  of  the  PWAFC 
be  in  its  lower  right  corner  (looking  from  cathode  side). 
Here,  the  gas  pressure  is  highest  and  the  circulation  of  liquid 
electrolyte  through  the  matrix  is  low,  as  shown  in  Fig.  21. 

In  few  cases,  as  explained  in  the  introduction,  we  have 
observed,  after  post-mortem  analysis  of  the  PWAFC,  that  the 
cathode,  anode  and  matrix  were  melted  together.  The  melted 
zone  was  repeatedly  observed  in  the  lower  right  region  of  the 
matrix  and  electrode  assembly  (looking  from  the  cathode 
side).  If  this  picture  is  superposed  on  the  picture  of  stationary 
pressure  distribution  on  cathode  (Fig.  22),  one  can  see  that 
the  damaged  zone  corresponds  to  the  high  gas  pressure  on 
the  electrode  and  to  the  relatively  low  pressure  of  the 
electrolyte,  exactly  with  the  zones  in  Fig.  21,  delineated 
by  ellipse. 


The  explanation  that  the  observed  damages  are  the  con¬ 
sequence  of  bad  construction  of  gas  and  electrolyte  distri¬ 
butors  and  not  the  consequence  of  cracks  in  the  electrode 
catalyst  layer  is  even  more  evident  when  the  superposition  of 
the  damaged  zone  is  made  on  the  picture  of  the  whiteprint 
paper  coloring  reaction  propagation  measured  with  the 
highest  flow  rate  of  ammonia  (Fig.  23).  The  observed 
coincidence  leaves  no  doubts  about  the  usefulness  and 
powerfulness  of  the  CFD  modeling  in  the  design  of  the  fuel 
cell. 


Fig.  23.  Superposition  of  the  photocopied  damaged  zone  after  post¬ 
mortem  PWAFC  analysis  on  the  propagation  contours  of  the  whiteprint 
paper  coloring  reaction  at  highest  ammonia  flow  rate. 
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7.  Conclusions 

The  CFD  study  of  the  mass  and  heat  transfer  in  the 
PWAFC  enabled  us  to  improve  the  design  of  the  fuel  cell 
in  order  to  achieve  as  uniform  as  possible  the  reactant  gases 
lateral  distribution  on  the  electrodes  and  heat  generation. 
This  is  important  because  the  ratio  between  the  lateral 
electrode  dimension  (in  X-  or  T-axis  direction)  and  the 
thickness  of  the  matrix  and  electrode  assembly  is  over 
50.  The  results  have  undoubtedly  pointed  at  changes  to 
be  made  in  the  following  details  of  PWAFC  construction: 

•  The  results  of  numerical  simulation  as  well  as  the  results 
of  experiments  with  ammonia  flow  in  the  cathode  com¬ 
partment  have  revealed  that  the  existing  geometry  of  the 
distribution  channels  for  gas  and  electrolyte  in  the 
PWAFC  is  not  proper.  Due  to  mismatch  of  electrolyte 
and  gas  flows  the  excess  heat  is  generated  in  some  critical 
points  which  lie  around  the  gas  inlet  point,  even  when  the 
fuel  cell  works  at  low  power.  When  the  power  increases 
this  phenomenon  becomes  more  expressive.  The  damages 
appear  locally,  like  that  shown  in  Fig.  22,  which  lower  the 
efficiency  of  the  fuel  cell.  In  order  to  avoid  this,  the  step¬ 
like  profile  of  the  gas  inlet  channels  must  be  eliminated, 
although  this  means  the  introduction  of  an  additional  layer 
and  the  increase  of  the  fuel  cell  thickness. 

•  The  cross-section  of  the  horizontal  channels  must  be 
increased  in  order  to  diminish  the  pressure  drop  along 
them. 

•  The  inlet  and  outlet  of  gases  and  electrolyte  must  be 
provided  from  both  sides  of  the  horizontal  channels  in 
order  to  eliminate  the  diagonal  distribution  of  gas  and 
electrolyte  on  the  electrode  surface  and  within  the  matrix, 
respectively. 

•  The  tightness  of  the  fuel  cell  must  be  improved  in  order  to 
prevent  mixing  of  electrolyte  and  gas  flow  paths. 

•  The  liquid  cooling  system  must  be  introduced  into  the  fuel 
cell  stack  design.  In  order  to  keep  the  temperature  changes 
of  PWAFC  at  peak  power  conditions  within  the  limits  of 
7-9°C,  we  would  need  the  mass  flow  of  cooling  water 
equal  to  0.11  kg  min-1  per  cell.  At  normal  working 
conditions  of  PWAFC  (current  density  about  0.3  A  cm-2) 
the  mass  flow  of  water  needed  would  be  4.7  times 
lower. 
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